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Abstract 

We consider the optimal Skorokhod embedding problem (SEP) given full marginals 
over the time interval [0,1]. The problem is related to the study of extremal martingales 
associated with a peacock (“process increasing in convex order”, by Hirsch, Profeta, 
Roynette and Yor [16]). A general duality result is obtained by convergence techniques. 

We then study the case where the reward function depends on the maximum of the 
embedding process, which is the limit of the martingale transport problem studied 
in Henry-Labordere, Obloj, Spoida and Touzi [13]. Under technical conditions, some 
explicit characteristics of the solutions to the optimal SEP as well as to its dual problem 
are obtained. We also discuss the associated martingale inequality. 

Key words. Skorokhod embedding problem, peacocks, martingale inequality, mar¬ 
tingale transport problem, maximum of martingale given marginals 

1 Introduction 

For a given probability measure /x on M, centered and with finite first moment, the 
Skorokhod embedding problem (SEP) consists in finding a stopping time T for a Brow¬ 
nian motion W, such that Wr ~ /x and the stopped process (MVa-) is uniformly inte- 
grable. We consider here an extended version. Let (nt)te[ o,i] be a family of probability 
measures that are all centered, have finite first moments, and are non-decreasing in 
convex order, i.e. t i-> := / R <)>(x)gt((lx) is non-decreasing for every convex 

function cf> : M —> M. The extended Skorokhod embedding problem is to find a non¬ 
decreasing family of stopping times, (A) t G[o.i]> for a Brownian motion W, such that 
W Tt ~ /Xi,Vt € [0,1], and each stopped process (MVla-) is uniformly integrable. We 
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study here an optimal Skorokhod embedding problem which consists in maximizing a 
reward value among the class of all such extended embeddings. 

For such a family // = (nt)t£ [o,i]> it follows from Kellerer’s theorem (see e.g. Kellerer 
[24] or Hirsch and Roynette [17]) that there exists at least one (Markov) martingale 
whose one-dimensional marginal distributions coincide with //. Assume in addition that 
t i '? m is right-continuous, then any associated martingale admits a right-continuous 
modification. Moreover, Monroe [26] showed that any right-continuous martingale can 
be embedded into a Brownian motion with a non-decreasing family of stopping times. 
This implies that the collection of solutions to the extended Skorokhod embedding 
problem is non-empty. Furthermore, the above optimal SEP is thus related to the 
study of extremal martingales associated with peacocks. A peacock (or PCOC “Pro¬ 
cessus Croissant pour l’Ordre Convexe”) is a continuous time stochastic process whose 
one-dimensional marginal distributions are non-decreasing in convex order according 
to Hirsch, Profeta, Roynette and Yor [16]. Since Kellerer’s theorem ensures the exis¬ 
tence of martingales with given one-dimensional marginal distributions, the interesting 
subject is to construct these associated martingales; we refer to the book [16] and the 
references therein for various techniques. We also mention that when the marginal 
distributions are those of a Brownian motion, an associated martingale is also referred 
to as a fake Brownian motion; see e.g. [1, 10, 20, 29]. 

Our problem to find the extremal martingales associated with a given peacock 
is motivated by its application in finance. Specifically, given the prices of vanilla 
options for all strikes at a maturity, one can recover the marginal distribution of the 
underlying at this maturity (see e.g. Breeden and Litzenberger [5]). Taking into 
consideration all martingales fitting these marginal distributions, one then obtains 
model-independent bounds on arbitrage-free prices of exotic options. The problem 
was initially studied using the SEP approach by Hobson [18] and many others. This 
approach is based on the fact that any continuous martingale can be viewed as a time- 
changed Brownian motion; let us refer to the survey paper of Obloj [27] and Hobson 
[19]. More recently, it has also been studied using the so-called martingale transport 
approach introduced in Beiglbock, Henry-Labordere and Penkner [4] and Galichon, 
Henry-Labordere and Touzi [11], Since then, there has been an intensive development 
of the literature on martingale optimal transport and the connection with model- 
free hedging in finance. In the present context of full marginals constraint, Henry- 
Labordere, Tan and Touzi [14] considered reward functions satisfying the so-called 
martingale Spence-Mirrlees condition, and solved an example of martingale transport 
problem with quasi-explicit construction of the corresponding martingale peacock and 
the optimal semi-static hedging strategy. 

In this paper, we study extremal martingale peacocks using the optimal SEP ap¬ 
proach. First, taking the limit of a duality result for a general optimal SEP under 
finitely-many marginal constraints, established in Guo, Tan and Touzi [12] (extending 
a duality result in Beiglbock, Cox and Huesmann [3]), we obtain a general duality 
result for the optimal SEP under full marginal constraints. Thereafter, we study the 
case where the reward function depends on the realized maximum of the embedding 
process. For the problem with finitely-many marginal constraints, the optimal ern- 
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bedding is then given by the iterated Azema-Yor embedding proposed by Obloj and 
Spoida [28], which extends the embeddings of Azerna and Yor [2] and Brown, Hobson 
and Rogers [6]. The solution to the associated dual problem, as well as the optimal 
value, is studied in Henry-Labordere, Obloj, Spoida and Touzi [13]. By applying lim¬ 
iting arguments, we obtain some explicit characterization of the optimal value and the 
primal and dual optimizers for the corresponding optimal SEP under full marginal 
constraints. 

The rest of the paper is organized as follows. The main results are presented in 
Section 2: in Section 2.1 we formulate our optimal SEP given full marginals, in Section 
2.2 we provide the general duality result, in Section 2.3 we focus on a class of maximal 
reward functions for which we specify the value of the problem and give the explicit 
form of a dual optimizer, and in Section 2.4 we present an associated martingale 
inequality. In Section 3 we provide further discussion of our results and relate them 
to the finite-marginal SEP. Specifically, we show that our full marginal optimal SEP 
has the interpretation as the limit of certain optimal SEP and martingale transport 
problems under finitely many marginals. The proofs are completed in Section 4. 


Notation, (i) Let 12 := C(K+,M) denote the canonical space of all continuous paths oj 
on R + such that uiq = 0, B be the canonical process, Po be the Wiener measure under 
which I? is a standard Brownian motion, F° = denote the canonical filtration 

generated by B, and F a = (F°')t>o be the augmented filtration under Po- 

We equip 12 with the compact convergence topology (see e.g. Whitt [31] or Stroock 
and Varadhan [30]): 


P(u,u') 


= E 

n> 1 


1 sup 0 < f < n | cut - Uj't 
2 n 1 + sup 0 < t < n I ut-u't 


Vcc, uj G H. 


( 1 . 1 ) 


Then (12, p) is a Polish space (separable and complete metric space). 

(ii) Let V+ = V+([0,1], R + ) denote the space of all non-decreasing cadlag functions 
on [0,1] taking values in M + . Similarly, let = V]*"([0,1],M) denote the space of all 
non-decreasing caglad functions on [0,1] taking values in M. 

Further, we equip V+ and Vj + with the Levy metric: for all 9 , 9' £ V+, 

d(M') := inf {e > 0 : 6 t _ £ - e < 0{ < 0 t+e + e, Vi € [0.1] }, (1.2) 

where we extend the definition of 6 to [— e, 1 + e] by letting 9 S := 9q for s £ [—e, 0] and 
9 S := 9 1 for s £ [1,1 -|-e]. Then V+ and Vj*" are both Polish spaces. 

(iii) As in El Karoui and Tan [8, 9], we define an enlarged canonical space by fl := 
12 x V+, where the canonical process is denoted by B = (B,T). The canonical filtration 
on the enlarged canonical space is denoted by F = (Ft)t> o, where Ft is generated by 
(B s ) se [ 0i t] and all the sets {T r < s} for s £ [0, t] and r £ [0,1]. In particular, all the 
canonical variables (Tr)re[o,il are F-stopping times. 

We notice that the a-field F oo := Vt>o coincides with the Borel cr-held of the 
Polish space 12 (see Lemma A.l). 

For a set V of probability measures on 12, we say that a property holds V— quasi- 
surely (q.s.) if it holds P— a.s. for all ¥ £ V. 
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(iv) Let Cb denote the space of all bounded continuous functions from R to R, and by 
Ci the space of all functions / : R —» R such that ppjyj € Cb- 


2 Main results 


Throughout the paper, we are given a family of probability measures on R, p = 
(/ i t)te[o,i]> satisfying the following condition. 

Assumption 2.1. The family of marginal distributions, p = (pt)t£[o, ip satisfies: 


|x| Ht{dx) < oo and 


x pt{dx) =0, t G [0,1]. 


Furthermore, po = <5{o}? t nt is cadlag w.r.t. the weak convergence topology, and p 
is non-decreasing in convex ordering, i.e. for every convex function 4> : R —> R, 


Fs{4>) < hM) ■= 


4>(x)p t (dx) 


whenever s < t. 


2.1 The optimal SEP given full marginals 

First, let us formulate the optimal SEP given full marginals. Let 'P(P) denote the 
collection of all Borel probability measures on the canonical space P, and define 

V := jp G 'P(iT) : B is a F-Brownian motion and 

B.^Ti is uniformly integrable under p|. (2-1) 

For the given marginals /j = (//t)te[o,il, we then define 

V{n) ■-= {PeiP :B Tt Vt e [0,1]}. 

Lemma 2.2. Suppose that Assumption 2.1 holds true, then V(/i) is non-empty. 

Proof. Since the marginal distributions p = (pt)te[ o,il satisfy Assumption 2.1, it 
follows from Kellerer’s theorem (see e.g. Kellerer [24] or Hirsch and Roynette [17]) 
that there is a martingale M such that M t ~ p t for all t € [0,1]. Since t 4 fij is 
right-continuous, the martingale M can be chosen to be right-continuous. It follows 
from Theorem 11 in Monroe [26], that there is a Brownian motion W and a family 
of non-decreasing and right-continuous stopping times (Tt) te r 0 ,ii , such that (W T1 a-) is 
uniformly integrable and {W T .) has the same finite-dimensional distributions as (M.). 
In consequence, the probability induced by (' W.,t.) on P belongs to V(p). □ 

The main objective of the paper is to study the following optimal Skorokhod Em¬ 
bedding Problem (SEP) under full marginal constraints: 

P(p) = sup E P [T(R.,T.)], (2.2) 

P eV(ti) 
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where <h : Q —> R is a reward function which is assumed to be upper semicontinuous 
and bounded from above. 

The optimal SEP (2.2) under full marginal constraints is given in a weak formula¬ 
tion. We specify this next. For given marginals /./ = (nt)te [o,ib let a //-embedding be 
a term 

a = {n a ,J^,F a ,¥ a = (^t>o,(Wnt>oAV) s em), (2-3) 

such that in the filtered space (D“, J 7 ", P Q , F“), W a is a Brownian motion, T° is a non¬ 
decreasing cadlag family of stopping times, the stopped process A .) is uniformly 
integrable, and W^ a ~ p “ fit for every t € [0,1]. Denote by A(/x) the collection of all 
/i-embeddings a. It is clear that every //-embedding a € A(//) induces on the canonical 
space II a probability measures P € P(//). Further, every P € "P(//) together with the 
canonical space D and canonical process B, forms a stopping term in A(//). It follows 
that the set P(/i) is the collection of all probability measures P on D, induced by the 
embeddings a € Al(//). As a direct consequence, the optimal SEP (2.2) admits the 
following equivalent formulation: 

P(n) := sup E" [<I>(lF a , T a )]. (2.4) 

2.2 Duality for the full marginal SEP problem 

In order to introduce the dual problem, let L^ oc denote the space of all F-progressively 
measurable processes, H = (H t )t> o, defined on the enlarged canonical space P and 
such that 

Hgds < oo, P-a.s., for all f > 0 and P € V. 

For every H € L f oc and P £ V, the stochastic integral of H w.r.t. the canoni¬ 
cal process B under P, denoted by (H ■ B)., is well-defined. An adapted process 
M = (Mt)t >o defined on O is called a strong supermartingale under P, if M r is inte¬ 
grable for all F-stopping times r > 0, and for F-stopping times t\ < r 2 we have that 
F p [M T2 \J 7 Ti ] < M T1 . We then define T~L by 

Ti := £ hf oc : {H ■ B). is a strong supermartingale under every P £ pj. 

Next, let M([0,1]) denote the space of all finite signed measures on [0,1]. Note that 
M([0,1]) is a Polish space under the weak convergence topology. Further, let A denote 
the space of all A : R —> Af([0,1]) admitting the representation 

A (x,dt) = Ao(x, t)X(dt), 

for some finite positive measure A £ Af([0,1]) and some locally bounded measurable 
function Ao : R x [0,1] —)• R. For // = (//t)tg[o,i]) we define 

A(//) := |a € A :/z(|A|) := [ [ \X 0 (x,t)\nt{dx)X(dt) < 00 ), 

1 Jo Jr j 
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and 


n(\) := I A (x,dt)fit(dx) = f f \o(x,t)nt(dx)\(dt), VA € A(//). 

Jo Jo Jr 

With the notation A (B) := [ ( ] Ao(L>t s , s)A(ds), we finally set 


(2.5) 


V(fi) := {(A, if) € A(n) x n :X(B) + (H-B) Ti > $ (B.) , P-q.s.}. (2.6) 


The dual problem for the optimal SEP (2.2), under full marginal constraints, is then 
defined as follows: 


D(») 


inf n( A). 

(A ,H)eT>(n) 


(2.7) 


Our first main result is the following. 


Theorem 2.3. Let Assumption 2.1 hold true. Suppose in addition that : Ll —> R is 
upper semicontinuous, bounded from above and satisfies $>(ui,0) = $>(u>q iA ., 0) for all 
(' uj,9 ) £ LI. Then, there exists a solution P € V(p) to the problem P(n) in (2.2), and 
we have the duality 


E f 


$(5.,T.) 


P(T) = D(ti). 


We also introduce the weaker version of the dual problem: 


Doth) 


inf 

(X,H)ev 0 (L) 




( 2 . 8 ) 


with Vq(h) given by 

:= {(A, AT) € A(/i) XU :\(B) + {H-B) Ti >*(B.), P(/x)-q.s.}. 
As a consequence of Theorem 2.3, we have the following result. 

Corollary 2.4. Under the same conditions as in Theorem 2.3, it holds that 


P (M) = A)M = D(n). 

Proof. Let (A, H) £ T>o (//■). For any P £ V(n), taking expectation over the inequality 
in the definition of T>q{h), one obtains //(A) > E p [<L(i?., T.)]. Hence, fi{ A) > P(n), 
which yields the weak duality Dq(h) > P(ff). Since V(jT) C V, it follows that Dq(h) < 
D(p). In consequence, the result follows from Theorem 2.3. □ 


2.3 Maximum maximum given full marginals 

In this subsection, we restrict to the case where 

4>(c jv, 6) = (,b{ojg j), with Wj := max u s , t > 0, (2-9) 
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for some bounded, non-decreasing and upper semi-continuous (or equivalently cadlag) 
function <f> : R + —> M. According to Theorem 2.3 and Lemma 4.3 below, we have the 
duality P(p) = D(p). Our main concern in this part is to compute this optimal value 
and, in turn, find and characterize the solution to the dual problem (2.8). 

First, we introduce some further conditions on the marginals /x. To this end, let 
c(t,x) := f R (y - x) + nt(dy) for every (t,x) G [0,1] X R. 

Assumption 2.5. (i) The function c is differentiable in t and the derivative function 
dtc is continuous, i.e. dtc(t,x) G C([0,1] x R). 

(ii) There exists a sequence of discrete time grids (7r n ) n >i with ir n = (0 = Vq < tf < 
■ ■ ■ < t™ = 1), such that |7r n | — > 0, and, for all n > 1, the family of finite marginals 
satisfies Assumption © in [28]. 


We introduce a minimization problem for every fixed m > 0. With the convention 
that § = 0 and g = oo for c > 0, let 


C(m) 


inf 

(GY]: C < m 


c( 0 ,Co) + f 1 d t c(s,Cs) dc 
m - Co Jo rn- ( s 


( 2 . 10 ) 


Our hrst result is on the value of the optimal SEP (2.2). 


Theorem 2.6. Let be given by (2.9) for some bounded, non-decreasing and cadlag 
function f. Suppose that Assumptions 2.1 and 2.5 (i) hold true. Then 

roo 

P{n) = D(n) < (j>(0) + / C(m)d(j)(m). (2-11) 

Jo 

Suppose in addition that Assumption 2.5 (ii) holds true, then equality holds in (2.11). 


Our second result is on the existence and characterization of a specific dual opti¬ 
mizer. To this end, for m > 0 and £ G V) + such that £(1) < m, let the associated 
functions A i' m and A^’ m be given by 


and 


x d 


(x,t) 


A i’ m (x,t) 


m — x 
{m - (t ) 2 




1 / (x - Ct) + 

A£t \ m ~ Ct 


(x-Ct+) + \ 
m-(t+ ) 


+ 


- Ci) + 

m- Ci 


1 P=i}> 


where := Ct+ ~ C t, Dm ■= {t € [0,1) : A( t > 0} and D r rn := [0,1) \ D m . We then 
define 

A <’ m (x,dt) := (Xi’ m (x,t) + X ( d ' m (x,t))d( t . (2.12) 

It is clear that A < ’ ,m G A. Next, we define the dynamic term. Let T m (ui) = inf{f > 0 : 
wg( t ) > m}. Further, let 0” 1 : M + [0,1] the right continuous inverse function ofsi-> 
9 S , given by Of 1 := sup{r G [0,1] : 0 r < s}. We note that 6 , F 1 (o)) is ^-measurable for 
fixed s and, thus, it is cadlag and F-adapted and therefore F-progressively measurable. 
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With I : R + —> M + and I + : M + —»• R + given, respectively, by I (s) := 6(6 1 (s)—) 
and I + (s) : = 6(6~ 1 (s)), we let 




m -Ce- 1 (Tm) r n~Ce~ 1 (s) 


(2.13) 


Finally, let £ : [0,1) x M + — > M such that, for all rn > 0, Q m € V+ and Cff 1 < rn, where 
C m = ((■,m ). Assuming that / 0 °° < oo, we then define 


/*oo /*oo 

A ^(x,dt) := / A 1, ’ m (x,dt)dcf>(m) and iff := / Rf ’ m d<f>(m). (2-14) 

Jo " Jo 

The construction of the dual optimizer below is based on the existence of a solution 
to the minimization problem (2.10). 


Lemma 2.7. Let Assumptions 2.1 and 2.5 (i) hold true. Then there exists a measur¬ 
able function ( : [0,1) x R + —>■ M such that, for all m > 0, ( m € V z + is a solution to 
( 2 . 10 ). 


Theorem 2.8. Suppose that <f> : M+ —>• M is non-decreasing and that Assumptions 2.1 
and 2.5 hold true. Let C : [0,1) x M + — >■ R be a measurable function such that, for all 
m> 0, C n e V+ is a solution to (2.10) , and 


dcj)(m ) 


I o (m — CD 2 


< oo. 


(2.15) 


Then, (A, H) := (A^, H^) G A(/x) x TL. Suppose in addition that <j) : M + —> M is bounded 
and continuous and that, for all t G [0,1], 


Ht is atomless; Q n and its inverse are both continuous in m. (2-16) 

Then (A ,H) is a dual optimizer for the problem Dq({jl) in (2.8). That is, with <f> given 
in (2.9), it holds that 

n(\)=D 0 (n) and A (B) + (H ■ B) Ti > $(R.,T.), V(/i)-q.s. (2.17) 

Remark 2.9. The condition (2.16) is needed to argue the convergence to (A ,H), in 
an appropriate sense, of the corresponding dual optimizers for the finite marginals case 
(see Lemma 4-4)- As seen from the proof, if Q m can be represented as a countable sum, 

i.e. 

OO 

C = £CW + ,]M, (2.i8) 

k=0 

for some (C*P)fc> 0 ; then (A ,H) is a dual optimizer even though condition (2.16) fails. 



2.4 An associated martingale inequality 

In this section we establish a closely related martingale inequality. We stress that this 
result does not require Assumption 2.5. 


Proposition 2.10. Let M be a right continuous martingale, <f> : R + —» R non¬ 
decreasing and cadlag, and £ : [0,1] x R + —>• R such that, for m > 0, £ m € V ; + 
and Cf 1 < m. Assume further that for each m > 0 ; 

|t € [0,1) : x i — » P [M t < x] is discontinuous at x = £f l j 

is a d£f l,c -null set, where £™’ c is the continuous part oftt-tQ 1 . Then, with Mf := 
maxo< s <(M S) it holds that 


e [mi)] 


< m + 




X m (M, dt) 


d(f>(m ), 


where X m (x,dt) := ,m (x, dt) (cf. (2.12)) so that 


E 


Uo 


X m (M, dt) 

e [(Ah - cr 


+ 


m - cr 


P [1AI t > CD (m - CD - E [(Mt - CD + ; 


+ 


E 


(m - Cj 
E [(M t - C 


m\ 2 
t ) 


d£t 


m\ 

t ) 


m — cr 


E[(Mt-cp + ; 

m — cr 


We conclude this section with a remark on an alternative version of the above 
martingale inequality. 


Remark 2.11. Suppose that M is a cadlag martingale such that the function CM(t , x) := 
E[(Mt — x) + ] is C 1 in t. Further, let £> : R + R non-decreasing and cadlag, and 
C : [0,1] x R + — > R such that, for m > 0, £ m £ ¥ z + and Cf 1 < m. Then, 


E 


mi 


< m + 



/ e[(Mq-CqD + ] 
V rn- Co m 



dtCMjtXf 1 ) 

m — cr 



d<f(m). 


Indeed, due to Monroe [26], there is P £ V(n) such that 


e [mi)\ 


E^ 


max Bt, 

0<t<l 


< E p [f(B*T 1 )]^ 


(2.19) 


where the inequality follows as <f is non-decreasing and maxo<t<i Bx t < Bff . The 
above inequality is therefore an immediate consequence of Theorem 2.6. 


3 Further discussion 

In this section, we provide some discussion on the relation between the optimal SEP 
and the martingale transport problem, and specify how the optimal SEP given full 
marginals can be considered as the limit of the approximating problem defined by a 
finite subset of marginals. We also discuss a numerical scheme for the problem C(m) 
in (2.10). 
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3.1 The optimal SEP given finitely many marginals 

First, we consider the optimal SEP given finitely many marginals and recall some 
results established in previous works. To this end, for n > 1, let 7r n = {ip, • • • , t r r \ } be 
a discrete time grid on [0,1] such that 0 = 2 q < if < ■ ■ ■ < 2” = 1. Then, let 

Vjj-i) ■■= |F G V : B Tf „ ~ P fj,q, Vk = 1, • • • ,raj. 

The set V n (ji) consists of all Skorokhod embeddings of the n marginals {lAvi)k=\,—n- 
Let <E> n : 12 x (M + ) n —> R be a reward function. The associated optimal SEP is 
formulated as 

Pn{n) := sup E p 
PePnOd 

3.1.1 The duality result 

In Guo, Tan and Touzi [12], a duality result is established for the optimal SEP (3.1). 
Let us define 

n 

D n (n) := inf j ^ //t«(A k ) ■ (Ai, • • • , A n , H) € (Ci) n x H such that 

k=\ 

U 

^A k{B Tt n) + (H ■ B) Ti > ® n (B. ATl ,T t n,---,T t n), P-q.s.|.(3.2) 

k=l 

One of the main results in [12] is the following duality result, which is a cornerstone 
in our proof of Theorem 2.3. 

Proposition 3.1. Suppose that Assumption 2.1 holds true and that <!?„, is upper semi- 
continuous and bounded from above. Then, P n (p 0 = B) n (ji) and the supremum of the 
problem P n (p) in (3.1) is attained. 

3.1.2 Optimal SEP and martingale transport problem 

One of the main motivations for studying the optimal Skorokhod embedding problem is 
the fact that any continuous local martingale can be seen as a time changed Brownian 
motion. It is therefore natural to relate the optimal SEP to the martingale transport 
(MT) problem. 

Let 12 := C([0,1],M) denote the canonical space of all continuous paths on [0,1], 
with canonical process X and canonical filtration F = (T))o<«i. Let A4 denote 
the collection of all martingale measures on 12, i.e. the probability measures P on 
(12, J 7 ) under which X is a martingale. We recall that there is some non-decreasing 
F-progressively measurable process (. X) which coincides with the quadratic variation 
of X under every martingale measure P G A4 (see e.g. Karandikar [23]). Let 

(X)f 1 := inf(2 : (X) t > a}. 

Then, under every P € M, the process (AT^-i ) s >o is a Brownian motion, and for every 
t > 0, (X)t is a stopping time w.r.t. the filtration (-P)x} _1 Let /j = (nt)o<t<i 


&JB. 


•ATi, 


It" 


(3.1) 
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be the given family of marginals satisfying Assumption 2.1. For n > 1 and a discrete 
time grid n n : 0 = < ■ ■ ■ < = 1, we denote 

M n (p) ■= {P € M :X t nr^n t n,k = l,---,n). 

For a reward function £ : fl —>• M, we then define the MT problem 

Pn{») ■■= sup E%(X.)1. (3.3) 

PG Mn(fl) 

The problem has a natural interpretation as a model-independent bound on arbitrage- 
free prices of the exotic option £(X.). In order to introduce the corresponding dual 
formulation, let TL denote the collection of all F- progressively measurable processes 
H : [0, 1] x — > M such that f H s dX s is a super-martingale under every F G M. c . 
Then, 

n 

D n (h) := inf | : (A, H) G (C\) n x TL such that 

fc=i 

n 

k (X t n) + (H-X) l > £(X), M-q.s.}. 

k =1 

The above dual problem can be interpreted as the minimal robust super-hedging cost 
of the exotic option, in the quasi-sure sense, using static strategies A and dynamic 
strategies H. 

Via the time change argument, the above MOT problem and its dual version are 
related, respectively, to the optimal SEP P n {n) and the associated dual formulation 
D n ([P). The following result is given in [12], It allows us to relate the limit of the 
above problems to our full marginal SEP; see Section 3.2 below. 

Proposition 3.2. Suppose that Assumption 2.1 holds true, and that the payoff func¬ 
tion £ : Cl —> M is given by 

ax.) = *n(x {x) - i A1 ,(X) t n,... ,(X) t n), (3.4) 

for some <f> n which is upper semi-continuous and bounded from above. Then, we have 
Puff) = Pn(h) = Dn(h) = D n (n). 

3.1.3 The iterated Azema-Yor embedding 

An example of payoff function £ satisfying the conditions in Proposition 3.2, is given by 
£(X.) := 4>{Xi) with X± := maxo<t<i X t and f : M + —» R a non-decreasing, bounded 
and cadlag function. This corresponds to the function defined in (2.9), for which 
the optimal SEP, given finitely many marginals, is solved by Henry-Labordere, Obloj, 
Spoida and Touzi [13], and Obloj and Spoida [28]. 

To solve this problem, a first technical step is to establish the following path-wise 
inequality (Proposition 4.1 in [13]). 
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Proposition 3.3. Let x be a cadlag path on [0,1] and denote x* := maxo< s <jx s . 
Then, for every m > xq and (1 < ■ ■ ■ < ( n < m: 


1 


( x t* n >rn} 


< 



( x q - Ci) + 
m - Ci 


+ l{x 



( ( x i; ~ 0+l) + . X *i+l x tA 

^ V m- 0+1 + {m - x v Ci+1 - Xti} m - 0+1 J ' 


(3.5) 


As argued in [13], the above inequality implies that also the following inequality 
holds: 

n rti 

l{x ( * n >m} < / H t’ m (x)dx t , (3-6) 

i=1 Jt i-\ 

with T m (x) := inf{t > 0 : x t > m} and 


A C,m 


(x) 



(x ~ Q) + _ (x - p+i)+ 
m-Ci m — O+i {l<n}l 
1 (^- 1 , t\ (T m fx)) + l[o,q_i] (-^m( x )) 1| X( i>0 } 

m - 0 


X' € M, 


^ ^ [^i—1 j 0) • 


Indeed, if x is continuous at T m (x), then the two inequalities coincide. If x has a jump 
at T m (x), then the first component of the dynamic term in (3.6) strictly dominates 
the corresponding term in (3.5). 

Intuitively, the l.h.s. of (3.5) can be interpreted as the payoff of a specific exotic 
option. It serves as the basic ingredient for more general exotic payoffs since any non¬ 
decreasing function 0 admits the representation p(x) := 0(0) + f Q l l x > m d(f)(m). The 
r.h.s. of (3.5) can be interpreted as a model-independent super-replicating semi-static 
strategy, the cost of which can be computed explicitly. 

Minimizing the super-hedging cost yields the following optimization problem: 


C n (m) 



C-i ( 0 + 1 ) 

m - 0+1 



(3.7) 


where Cfc(x) := /0 (y — x) + fit k {dy) . It is argued in [13] that the minimization problem 
(3.7) admits at least one solution (■ Ck( m ))i<k<n■ An immediate consequence is that 

poo 

D n (fi ) = D n {n) < 0(0) + / C n (m) d(f>(m). (3.8) 

Jo 

Under further conditions (Assumption © in [28]), Obloj and Spoida [28] provide an 
iterative way to solve (3.7), and to obtain a family of continuous functions ( fk)i<k<n 
satisfying (k(m) = max/j<j< n O(m), Vm > 0. Using the family of functions ( fk)i<k<n , 
they further dehne a family of iterated Azema-Yor embedding stopping times, given 
by tq := 0 and 




inf{t > r fc _i : B t < 4(5*)}, if B Tk _ x > 4 (B 
T k -i, else. 


(3.9) 
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The stopping times, {rk)k=i,- ,m embed the marginals {nt™)k=i,— ,n- Moreover, it is 
proven in [13] that the embedding satisfies 

POO 

E[0(W?J] = 0(0)+/ C n (m)d<f>(m). 

J o 

In consequence, under Assumption © in [28], it holds that 

POO 

P n {p) = P n {p) = P>n{p) = D n {n) = 0(0)+ / C n (m)d(j)(m). (3.10) 

Jo 

We notice that the discrete process {W Tk )i<k<n resulting from this construction is in 
general not a Markov chain. 

Remark 3.4. In [13], the result (3.10) is formulated for the continuous martingale 
problem as defined in (3.3). However, it can be easily deduced that the solution is 
optimal also for the corresponding cadlag martingale problem. Specifically, let Cl d de¬ 
note the space of all cadlag functions on [0,1], X the canonical space with canonical 
filtration F d , and M d the space of all martingale measures. Define 

P d {p) := sup E # [0(Xn], with M d (p) := {p € M d : X tk /^,Vfe}. 

It is clear that P n {p) < P d {p) since every continuous martingale is a cadlag martin¬ 
gale. Further, by Monroe's [26] result, every cadlag martingale can be represented as 
a time changed Brownian motion. Since maxo<t<i u$ t < Ug and 0 is non-decreasing, 
it follows that P d (p) < P n (/i). Therefore, according to (3.10), for the payoff f(X.) := 
4>(X[) with 0 : M —> R non-decreasing, 

PM = pM- 

3.2 The optimal SEP given full marginals 

Our optimal SEP (2.2) given full marginals is obtained as the limit of the problem 
given finitely many marginals; see the proof of Theorem 2.3. We provide here fur¬ 
ther discussion of the convergence of various optimal values and the corresponding 
optimizers. 

3.2.1 The limit of MT problem given finitely many marginals 

Our main motivation for studying the optimal SEP is the MT problem, which has 
a natural interpretation and applications in finance. For the case of finitely many 
marginal constraints, and for certain payoffs, the optimal SEP P n (p) in (3.1) is equiv¬ 
alent to the MT problem P n {p) in (3.3) (cf. Proposition 3.2). 

When the number of marginals turns to infinity, the question is whether the MT 
problem (3.3) converges in some sense. Specifically, we are interested in the conver¬ 
gence of the optimal value and of the optimizer. The following convergence result is 
an immediate consequence of the proof of Theorem 2.3. 
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Proposition 3.5. Suppose that Assumption 2.1 holds true and let £ : fl —» R given by 

«*■) = 

for some upper semicontinuous and bounded $ : x R + —> R. Let P n (jj) defined w.r.t. 
f in (3.3) and P n {h) an d -P(/f) defined w.r.t. $. Then, we have the approximation 
result 


lim P n (l-i) = lim P n {p) = P(fjt). 

n^-oo n—>oo 

Further, the optimal transferences converge in sense of the convergence of Skorokhod 
embedding (i.e. the convergence of probability measures on Ft). 

Remark 3.6. Recall that Q is the canonical space of continuous functions on [0,1], 
we define 

P(p) := sup E^[£(X)], with M(p) := jlP G M. : X s /i t t G [0,1]|. 

As we can see below, the limit of the optimal (n marginal) continuous martingales may 
be a cadlag martingale. Thus the convergence of the MT problem P n (h) to the MT 
problem P{n) fails in general. This underpins the importance of the full marginal SEP 
as the correct way of specifying the limit of the continuous n-marginal pricing problem. 


3.2.2 The limit of the optimal martingale transference plan 

We discuss here a specific case where the limiting martingale can be explicitly char¬ 
acterized. Specifically, Madan and Yor [25] provide, under certain assumptions, a 
characterization of the continuous time martingale obtained from the Azema-Yor em¬ 
bedding. Let bt(x ) be the barycenter function of p t , defined by 


b t {x) 


4,oo) yht(dy) ^ 

lk([x, oo)) x<r " 


+ xl X >r M ■ 


Suppose that the family (nt)te[ o,il satisfies the so-called property of increasing mean 
residual value: 


t !-)■ bt(x) is non-decreasing for every x. (3.11) 

For any discrete time grid 7r n : 0 = tp < t” < • • • < t” = 1, and under the additional 
Assumption © in [28], it turns out that the boundary functions ( Ck)i<k<n are then 
given by fk : = b -1 , and that the iterated Azema-Yor embedding coincides with the 
Azema-Yor embedding: 


T t := inf{s > 0 : B s > Ct(B* s )}. 

Notice that the iterated Azema-Yor embedding induces a continuous martingale, which 
is the optimal martingale transference given finitely many marginals. It follows that 
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under condition (3.11), this optimal martingale transference plan converges to the 
left-continuous right-limit martingale M = (Mt)te[o,ii, given by 

M t := B n . 

In [25], the authors prove that M is in fact a Markov process and provide its generator 
in explicit form. 

To conclude, we highlight that the Markov process M defined above is a left- 
continuous process, a right-continuous modification gives the same generator. It is 
easily verified that 

Mf < B* T1 and F[Ml < B * T J > 0. 

In consequence, M provides no solution to the MT given full marginals; see also 
Remark 3.6. 


3.2.3 The limit of the pathwise inequality 


The proof of Theorem 2.8 is based on applying limiting arguments to the path-wise 
inequality (3.5) (cf. Section 4.4). By use of a similar argument, we might obtain an 
almost sure inequality for cadlag martingales. 


Proposition 3.7. Let M a right-continuous martingale, f : [0,1) x R + —> R such that 
Q m G ¥/■ and Q m < m, and <f> : R + M + bounded, continuous and non-decreasing. 
Suppose either i) that (, (j) and the marginals of M satisfy the conditions of Theorem 
2.8; or ii) that f admits the representation (2.18). Then, with A^ given in (2.14), M 
satisfies the following inequality: 


<j)(MD < [ \<(ds,M s ) + 

Jo 


0 J 0 


#(m) 

{m<M*_ ;C r<M t - } m - ‘ ’ 


a.s. 


(3.12) 


The difference between the r.h.s. of (3.12) and (2.17) appear in the dynamic terms 
(cf. (2.13)). Specifically, for the martingale formulation, the counterpart of the first 
dynamic term in (2.17) is always negative and thus vanishes from the inequality. This 
is related to the fact that the limit of the first dynamic component in (3.6) is zero. For 
continuous martingales, the two inequalities coincide. 


3.3 The resolution of C{m) 


Finally, we would like to discuss the resolution of the problem C(m ) in (2.10), since 
the main results in Theorems 2.6 and 2.8 rely on its solution £. 

First, it is clear that we can decompose the minimization problem C(m ) as follows: 


C{m) = inf 


x<m 


c(0, x) 


m — x 


+ u(0, x) > ; u(0, x) := inf 


d t c(s, Cs 
m-Cs 


-ds. 


t£Vf,( o =xJ0 

The problem to compute v(0,x) is a standard singular deterministic control prob¬ 
lem. When the function d t c(s,x) is continuous, it therefore follows by standard argu¬ 
ments (see e.g. [7]) that v can be characterized as a viscosity solution to the PDE 

d t c(t, x) 


max < — d x v(t, x), — dfv(t,x) — 


m — x 


= o, 


(3.13) 
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equipped with the terminal condition x(l,x) = 0, for all x < nn. 

We now propose a numerical scheme for the problem C(m). To this end, for a 
given partition 7r n = {tf ,..., t f f }, with 0 = tg < ... < = 1, let ¥” the subset of ¥ ; + 

for which £ is constant on i = 1, ...,n. Further, let 

„.n/n ^ f ' $ C ( S >&) . 


u n (0,x) := inf [ 

Cevp, i 0 


m-Cs 


C,o=x, Ci < m 


For a sequence of partitions such that |7r n | —> 0, it follows that v n (0,x) —> u(0, x); cf. 
the proof of Lemma 4.2 below. On the other hand, 


x<(i<...<C„<m 


ft ’ 


with Ac(i”, C) := c(t”,£) — c(t”_ 1 ,C)- In consequence, u n (0,x) = u n (tg,x), where 
v n (t^,x), k = 0, ...,n, is iteratively defined by 

f «"(**, ®) = info < y<m - x {v n (tl +v x + y) + , k<n-l, 

\ v n (t%,x) = 0. 

This yields a scheme for explicit calculation of u n (0, x) as an approximation of n(0, x). 


4 Proofs 

4.1 Technical lemmas 

Lemma 4.1. Let <J> : hi —> R be upper semicontinuous, bounded and such that 4>(w, 6) = 
4>(w. A 6»ii 0), for all (uj,0) G Q. Further, let (7r n ) n >i be a sequence of discrete time grids 
with \w n \ —> 0. Then, there exists a sequence ($ n )n>i of upper semicontinuous and 
bounded functions <f> n : fl x (M + ) n —>• M, such that 

4>„(w. a ,o 1 ,Oq,<..,9 t n)\$(uj. A o 1 ,0.), asn-> oo, V(fl,w) £ H. (4.1) 

Proof. Let 4> fc : O x ¥+ —» M, k G N, such that <F fc (u;,0) = ^ k (oj. A g 1 ,0), <F fc is 
bounded and Lipschitz and \ <h. Then, let 4> fc,m : Q x ¥+ —> R, m G N, be given 
by 

$ k ' m (uj,0) := <L n (c o,0 m ), 

where 0 m := 0 ti , s G [^, ^), i = 1, for 0 G V+. Note that since d(0 , 0 m ) < ^ 
(cf. (1.2)), we have 

$ fc ’ m (u;,60 -$ fc (cu,(9) = $ fc (w, @ m ) - $ k (u, 0) < —, 

m 

with Lk the Lipschitz constant associated with <h fc . In consequence, 

4> fc ’ m (u;, 0) := 4> fc ’ m (w, 0) + — \ $ fc (cn, 0), m -> oo. 

m 

Hence, we may choose such that & k,mk (uj,0) \ 4>(cn,0), k — > oo. In consequence, 
defining <3? n := with k n = max{/c G N : < n}, we have that 4> n , n G N, 

satisfy (4.1) and we conclude. □ 
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Proof of Lemma 2.7. We follow the argument at the beginning of Section 3 of [13]. 
Let 


®m(C) 


c(0,Co) + f 1 fyc(sXsl ds 
m - Co Jo rn- C, s 


(4.2) 


We first consider a constant function ( z = z for some constant z < m. By direct 
computation, it is easy to see that 

'MC) = > C(m). 

m — z 

Note that since z i->- c(l, z) is convex and is the slope of the tangent to z c(l, z) 
intersecting the x-axis in rn, it follows that C(m ) < 1. 

On the other hand, since dtc(s, z) > 0, we have 

^m(C) 'i ~1 as Co ~> —oo. 
m - Co 

For the minimization problem C{m) in (2.10), it is therefore enough to consider the 
space ¥+([0, 1)> [K,m)) for some constant K £ (—oo, m), i.e. 


C\m) = inf M(C)- 

CGV+([ 0,1 ),[K,m)) 

Notice that ( i->- \H m (C) is continuous and ¥ z + ([0,1), [K, rn)) is compact under the 
Levy metric. It follows that, for every m > 0, there exists at least one solution in ¥+ 
to (2.10). 

To conclude, it is enough to use a measurable selection argument to choose a 
measurable function Q. □ 


Lemma 4.2. Recall that C n (m ) is defined by (3.7). Suppose that the function c is 
differentiable in t and that the derivative function dtc is continuous. Then, for every 
rn > 0, we have 

lim C n (m ) = C(m). 

n—>• oo 

Proof. Let Yf the subset of ¥ ; + for which £ is constant on (tj_i, tf\, i = 1,..., n — 1, 
and on (t n -\,t n ). For n fixed and ( £ Yf, let ((t n ) '■= C(tfi)- Notice that for every 
C€¥f, 

4KC) = + V [ U dtc{s ^ u) ds 

m - Cto Jti-1 m - c U 

c{Q,Ct 0 ) + y Cc(M)_MMJA 

™ ~Ct 0 ~[\ m - C u m -Cu ) 

_ v-^ ( Cu) _ c {U,Ct i+ J \ 

-C<> '"-Co-, i<n )' 
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Next, for any £ € V+, by direct truncation, we can easily obtain a sequence £ n 
such that £ n € V” and £ n —> £ under the Levy metric. It follows that C n (m ) —>■ C(m ) 
as n —> oo. □ 


Lemma 4.3. TTie mapping from 17 to M given by: 

(cu,9) i—» w*(#(l)) := sup w(s), (4-4) 

o<s<e(i) 

is continuous with respect to the product topology on 17. 

Proof. Let (w n (-), #”(•)) € 17, n € N, converging in the product topology to 
(&(•), #(•)) € 17. Recall that C(M+,M) is equipped with the metric p defined in (1.1), 
which induces the topology of uniform convergence on compact subsets. Hence, 

lim sup |u; n (s) — w(s)| = 0, for all m > 0. (4.5) 

n_5>oo 0 <s<m 

Further, convergence in the Levy metric is equivalent to point-wise convergence at 
each point of continuity. Due to the right-continuity of elements in V+([0,1],M+) and 
the extended definition of the Levy metric (cf. (1.2)), it follows that 9 n ( 1) converges 
to 9(1). Note that 


sup oj n (s) — sup u>(s) 

= sup w n (s) 

— sup Ld(s) 

O<s<0 n (l) 0<s<8(l) 

O<s<0 n (l) 

0<s<d n (l) 


+ sup D(s) 

— sup w(s) . 


O<s<0"(l) 

O<s<0(l) 


The first term is dominated by sup 0<s<e n( 1 ) |w n (s) — D(s)| which tends to zero as n 
tends to infinity due to (4.5). Since £>(•) is a continuous path, also the second term 
tends to zero. Hence, the mapping in (4.4) is continuous and we conclude. □ 

Lemma 4.4. Let tl d := Z7([0,1],M) be the space of all cadlag paths on [0,1] with 
canonical process X, and M d be the space of all martingale measures on & d . We 
define 

M d (p) := <jP€A4 d VtG [0,1]}. 

Further, let £ : [0,1) — > (—oo, m) be a non-decreasing caglad path on [0,1) and n n : 0 = 
7q < • • • < 7"; = 1 be a sequence of discrete time grids such that \-ir n \ —>• 0 as n ^ oo. 
Let £ c the continuous part of £ and let defined in (2.12). Then, 
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i) if fit is atomless, t G [0,1], we have M. d (/a)— q.s., 


n 


E 


( ( Xt l ~ ^) + 

V m - k 


( x q ~ Cw k+1 ) + 
nn — Ct n 

^ L k+1 


'-{k<n} 



ii) If Cs = J2T=o Ckl(t k ,t. k+1 ]( s )> then the convergence in i) holds path-wise for all 
x G fl d . The integral with respect to df c is then identically zero. 

Proof. It follows from the definition of X^ m , that in order to prove i), it is sufficient 
to show that, q.s., 


n— 1 


E 


( ~ ^+i) + 

V m - Ctjf +1 


-Cig) + \ 
m - Ct» ) 


X,- 


(Xt - Ct+) + 

m - (t+ 


(4.6) 


(X t ~Ct) + \ 

m - C, t ) 


Observe that for each path x6D([0,l],M), the discrete sum in (4.6) might be written 
as /g 1 f n {t ; C)dft, where 


fn(t;C) 


y 1 *e(«g,«g +1 ] ( ( x ^ ~C^ +1 ) + _ (x t n ~Cig) + \ 


Denote by D £ C (0,1) the subset of all discontinuous points of (. First, suppose that 
assumption i) holds. Then, for t 0 n {t : x t = Ct}j the /„(•; £) converges point-wise 
to /(■;<), with 


mo 


x-t—m i . 

(■ m-Ct) 2 ■ Lx ‘>0> 

1 ( (Xt-Ct+)+ _ (x«-Ct)+ \ 

c t +—Ct v m ~c*+ m ~c* y ’ 


t g 

i g D c . 


(4.7) 


On the other hand, by use of Fubini’s theorem and assumption i), we obtain that for 
P € M d (p), 


E p 




P[x t = Ct]dC t c 



Tt({Ct})d(f 


0. 


That is to say, 1 {x t =£t}d(t = 0, Xi d (fi)— q.s. Since, for all e > 0, £ —;> is 

Lipschitz on (—oo, m — e], there is K > 0, such that f n (t) < K,t G [0,1], n > 0. Hence, 
by use of dominated convergence we obtain /)' f n (t ; C)<7Ct —>■ [ ( ] f(t ; ()dQ, M d (p)—q.s., 
which implies (4.6). 

Next, suppose assumption ii) holds. Then, for t G [0,1], the function / n (t;£) 
converges point-wise to /°(t;£), where 


/°(t;C) 


o, 

1 / (xt-Ct+)+ _ (x«-Ct)+ \ 

c t +-c* v m_ c*+ m -Ct / ’ 


t G D'l 
t G T> c . 


(4.8) 


By use of the same arguments as in the case i), we may then apply the dominated 
convergence theorem pathwise and we easily conclude. □ 
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4.2 Proof of Theorem 2.3 


We first argue that the optimal SEP given finitely many marginals defined by (3.1), 
may be reformulated similarly to that in (2.3) and (2.4). Concretely, for a given 
discrete time grid 7r n : 0 = < ■ ■ ■ < = 1, we call a (/U, 7r n )-embedding a term 


a = ^f.PfP^F" = (7 r “)i> 0 , (W“)t>o, (T^)fc =li ... in ), (4.9) 

such that in the filtered space (14", I 7 ", P", F"), W a is a Brownian motion, T" < • • • < 
T“ are all stopping times, the stopped process (Wfk A .) is uniformly integrable, and 
~ p “ iM;n for each k = l,-- - , n. Let An{p) denote the collection of all (/x,7r n )- 
embeddings a. Then it is clear that every term in A n (/i) induces on the canonical 
space 14 a probability measure in V n (n), and every probability measure P £ V n (p\) 
together with the space (14,1 7 , F) forms a (^,7r n )-term in A n {p)- And hence, for 
<h n : 14 x (M+) n —> R a given reward function, we have that 


Pn(h') 


sup E p “ 
aeAi(At) 


*n(w a ,T 



(4.10) 


Before proving Theorem 2.3, we present a Lemma. Its proof is partly adapted from 
the proof of Theorem 11 in Monroe [26] and that of Theorem 3.10 in Jakubowski [22], 

Lemma 4.5. Let a n £ A n (n), n £ N, be a sequence of terms of the form (4.9). Let 
P ra be the probability measure on Q induced by (lL mn ,T“ 71 ) in the probability space 
(jl° ln ,J 7an , P an ). Then, the sequence {P n } n >i is tight, and any limiting point P is in 
V{p). 

Proof, (i) We first claim that the sequence {P«} n >i is tight. Indeed, the projection 
measure P n | q on Q is the Wiener measure for every n > 1, and hence the sequence 
(Pn|fi)n>i is trivially tight. Next, since T“ n are all minimal stopping times in the sense 
of Monroe [26], it follows from Proposition 7 in [26] that 

Pn (Ei > A) = P°" (If" > A) < A -1 / 3 (E p “ n [\Xf n |] 2 + l) , VA > 0. 

Let A\ be the set of functions in V p ([0,1], R + ) which are bounded by A > 0 and 7r the 
projection of 14 onto Vf ([0,1],M + ), it follows that 

P n (7T- 1 (A A )) = P n m<A) > A- 1 / 3 (( /Ul (|x|)) 2 + l). (4.11) 

Since A\, A > 0, are compact and the r.h.s. of (4.11) can be made arbitrarily small 
by an appropriate choice of A, the sequence of projection measures (P n | v +) n >i is also 
tight. In consequence, {Pn} n >i is tight. 

(ii) Let P be a limit point of (P n )n>i, by taking subsequences if necessary, we can 
assume that P n —> P. We now prove that B is a F-Brownian motion under the limit 
measure P. 

Since the measures P n are induced by (IT Q ",T Q ") under P“", we know that I? is a 
F-Brownian motion under each P n , n > 1. Let 4>s,0<e<4 — s, and 4> : 14 —> R be 
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a bounded continuous function which is T s+£ 
we have 


<f>(B.,T.){<p(B t )-<p(B s+e )- 




By taking the limit n —>• oo, it follows that 


E 


<f>(B.,T.)[<p(B t )-<p(B a+e )- 


-measurable, then for every <p € 


-(p"(B u )du 


' S+£ 


= 0. 


f \v"{B u )du 

>s+e A 


= 0. 


(4.12) 


According to Lemma A.l, the equality (4.12) holds true also for every bounded random 
variable 0 : fl —\ R that is ^-measurable. Let e —> 0, it follows that for every 
4> : 14 —> R bounded and ^-measurable, and every p G C%( R) that 


E* 


<j>(B.,T.) ( <p(B t ) — <p(B s ) - 


-<p"{B u )du 


= 0. 


And hence B is a F-Brownian motion under P. 

(iii) Now, we show that the process (-BtATi)t>o is uniformly integrable under P. For 
every e > 0, there is K e > 0 such that 



K e ) + m(dx) < e. 


Since |aj|l{|a;|> 2 A'} < 2(|a;| — K) + , it follows that 


E 


I / w 


< 2E " 


Bt i 


-K f 


+■ 


■ 2s, 


Then, for every bounded continuous function p : K —> M such that p(x) 
it follows by the dominated convergence theorem that 


E^ 


p(B t atJ 


lim E Pn 

n—¥ oo 


p(B t atJ 


< 2e, Vt > 


Vt > 0. 

< kl 1 (h|>2AV} 1 


o, 


which implies that (Bt/\Ti)t>o is uniformly integrable under P. 

(iv) Next, we prove that Bx t ~ p pt, t G [0,1]. We shall adapt the idea of proof of 
Theorem 3.10 in Jakubowski [22]. For every a n the solution of the optimal SEP (3.1), 
denote by m an the random measure on ([0,1],£>([0,1])) defined by 

m an {[ 0,t ] ,07) := 1+ * I Ll u y V ^[0,l], 

Since m n takes value in a compact space (the space of all positive measures on [0,1] 
with mass less than 1), we know that the sequence of distribution of m n under P“ n 
is tight. By taking subsequences and the Skorokhod representation theorem, we can 
assume that there is some probability space (fi*, J 7 *,P*) in which 


(' W an ,T? n ,m an ) -»• (W*,T*,m *), P * - a.s. 

Further, the map 1 1 — > E p * [m*([0,t])] from [0,1] to R is non-decreasing, and hence 
admits at most countable discontinuous points. It follows that there is some countable 
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set Qi C [0,1) such that E p * [m*({f})] = 0, for every t £ [0,1] \<Qh. Thus, for every 
t £ [0,1] \ Qi, we have P*-a.s., m an ([0,t]) —m*([0,t]), and hence T" 71 —>• T*. In 
particular, we have 


Wp n -> Wfy, P* - a.s. Vf £ [0,1] \ Qi. 

Besides, by Hirsch and Roynette [17] Lemma 4.1., there exists a countable set Q 2 C 
[0,1] such that t 1 —>■ /x* is continuous at any s £ [0,1] \ Q 2 . Then for every t £ 
[0,1] \ (Qi U Q 2 ), we have 

P o (B Tt )- 1 = P* o (Wj*) _1 = fi t . (4.13) 

By the right continuity of fit, it follows that (4.13) holds true for every t £ [0,1]. 

In summary, we have proven that in the filtered space (fi, J r 00 , P, F), B is a Brown¬ 
ian motion, T\ is a minimal stopping time and Rji ~ IM for every t £ [0,1]. We easily 
conclude. □ 


Proof of Theorem 2.3. By taking expectation over each side of the inequality 
defining in (2.6), for all P £ V(n), we easily obtain the weak duality P{n) < D(n). 
Let <h n , n £ N, the sequence of functions approximating as given in Lemma 4.1. 
Further, let P n {n) and D n (n) the primal and dual n-marginal problems defined w.r.t. 
in (3.1) and (3.2). Since <L n > 4>, it follows that 

D n (n) > D(n). (4.14) 


For each n £ N, let a n £ A n (ji) be the solution of the optimal SEP P n {n) in 
(4.10) defined with respect to <h n . Let P n be the probability measure on f! induced by 
(l'F Q ' a , T° n j in the probability space (fl an , P an Then, according to Lemma 4.5, 
the sequence {P n } n >i is tight, and P £ V(n), where P is a limiting point of {Pn} n >i- 
Note that by taking sub-sequences if necessary, we can assume that P n —>• P. By 
use of the monotone convergence theorem and the optimality of the P n for the problem 
(3.1), n £ N, it follows that 


P(H) > E P U(R.,T.) 
> lim ( lim 

n->oo 00 

= lim P n {n). 

n—>00 


= lim E p 

n—> 00 

$n(£.,T.) 


M R ’ T ) 


> lim ( lim E Pfe 

n—>00 V fc—>■ 00 




In consequence, since P n {n) > P(jJ-) for all n > 1, we have that 


lim P n (n) = P(n). (4.15) 

n—> 00 

Since the 4> n satisfy (4.1), we may apply the duality result for the optimal SEP with 
finitely many marginal constraints (see Proposition 3.1). Hence, it follows from (4.15) 
combined with (4.14) that P{n) > D(/i). Combined with the weak duality, this yields 
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P(fi ) = D(fi). As a by-product, we also obtain that P is an optimal embedding for the 
optimal SEP (2.2). This concludes the proof. □ 

Proof of Proposition 3.5. Given the form of £ : f2 — > M, Proposition 3.2 applies. 
Hence, P n {n) = Pn{Y)- Next, note that P n (Y) is of the form (3.1), for all n £ N. Let 
P ra be the optimal measure for P n ( / u). Then, according to Lemma 4.5, passing to a 
subsequence if necessary, P n —> P, with P £ V(n). It follows that 


P(fi) > E p $(B.,T.) 


> 


lim E Fn 

n—> oo 


§(B.,T.) 


Since P n (n) > P(/x), for n > 1, we easily conclude. 


lim P n (n). 

n—>oo 


□ 


4.3 Proof of Theorem 2.6 

By Lemma 4.3, the mapping 


(cj(-), 0(-)) i—> sup u(s), 

O<s<0(l) 


is continuous with respect to the product topology on ft. Hence, Theorem 2.3 applies 
and linx )WOO P„(^) = P(n) = D{(x). 

Next, for the optimal SEP with finitely many marginals, we have (see (3.8)) 


Pn{I- 1 ) 


D n (/J,) < ^(0) + 



C n {m)d<f)(m), 


where the equality holds under Assumption 2.5 (ii) (see (3.10)). 

Finally, it is enough to use Lemma 4.2 together with the monotone convergence 
theorem to deduce that 


/*oo 

lim / C n (m ) d(j)(m ) 

n^oo J 0 



C{m) d(j)(m). 


□ 


4.4 Proof of Theorem 2.8 


Due to assumption (2.15), the pair (A ,H) is well-defined and (A ,H) € A(/x) x H. 
According to (2.5) and (2.12), for m > 0 and ( £ V) + such that ( < m. the cost of 
is given by 


/r(A c ’ m ) = 

^ h 


c(0,Co) 
m - Co 


+ 


'c(t,Ct+) 

c(t,CtY 

r 1 d 

(c(un 

m-(t+ 

m - Ct_ 

Jo d( 

l m - C / 

1 d t c(s, (s 

- ds , 



ni - (s 





x—m 

>C1 (m—Q' z 

. Since ( m 

minimizes 




C=ff 


that /.i(A^ m,m ) = C(m). Integration w.r.t. d(j)(m ) and application of Theorem 2.6 and 
Corollary 2.4, yields //(A) = Dq(h). 
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Next, let 7r n : 0 = tfi < ■ ■ ■ < t™ = 1 be a sequence of discrete time grids such that 
7r n | —> 0, as n — > oo. According to Proposition 3.3, for {9,u) G fl, 


10 

If , 1 < V 


i u 9 (u) ~ CtT) + - CtT +1 ) + 

* * -*-{i<n} 


m - C™ 


m - CtT +1 


+£* 


m — 


( w , 


i=l 
n —1 

£* 

2=1 




<m<u* J 


^(u) 


m _ (~rn 
^ti 

u 0(t i+1 ) ~ Uo( ti ) 


{™<u* e(ti y <™ +1 < "«(**)> 


m - CT +1 


(4.16) 


Note that the l.h.s. does not depend on the partition. Hence, in order to verify that 
(A, H) satisfies (2.17), it suffices to show that, for all P € integrated w.r.t. d<f)(m ) 

the r.h.s. in (4.16) converges P-a.s. to A(Bt s ,cIs) + / Q Tl H s dB s . 

Application of Lemma 4.4 with X t = B t m and ( = ( m , yields that the static term 
in (4.16) converges to 


/V>» (B T ( t ), dt) , 

Jo 




In consequence, integrated w.r.t. d<j>{m ), the static term in (4.16) converges to A (B) = 
fo M B T s ,ds), P(fi)~ q.s. 

As for the first dynamic term in (4.16), the definition of r m yields: 


£ 

i=l 


m — 


-{p 






u 8(U) 


m-C 


m 

U 


( '° {L) 1 [r m ,/+(r m )]('S) 


m - C-i 


dB s , V — q.s. 


9 1 (r m ) 


Integration with respect to dcj), then gives the convergence of the corresponding terms. 

To prove the convergence of the second dynamic term in (4.16), we first integrate 
with respect to d<f). The integrated term may then be re-written as follows: 


n— 1 




d(j){m) 


{ u e(t i+1 ) -ueiu)) 


m - Q+i - w »(q)} m — Q 

( u »(.U)) Au e(.ti) d(j)(m ) r d di+i) 


H+l 


£ 


i —1 


r 0 


m ~ Cu +1 J9(ti) 


dB q 


^(i)” 1 d(j)(m) 


£ 

2=1 


m-C 


m 

U+i 




where denotes the inverse of Ct i+1 - We denote the integrand in (4.17) by H n . 
Note that the H n are predictable. Further, since ( satisfies (2.15) and is bounded, 
the H n are uniformly bounded. Due to assumption ii), we also have that H n —> H on 
D x (0, oo), with 


H s 




(s) \ i~M 


Au; 


■(*) 


d(j)(m ) 


(S ) 
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In consequence, for all P € V, we have that f^ 1> H™dB s —> f'^ 1) H s dB s in probability. 
Hence, convergence holds a.s. along a subsequence and we conclude. □ 


Proof of Proposition 3.7. Let m > 0 fixed and 4>(x) = l{ x > m y, the general case 
follows by integration with respect to dcj)(m). Recall that (3.5) holds for all cadlag 
paths x (in the proof of Theorem 2.8 we only made use of this result for continuous 
paths). For a given sequence of partitions n n : 0 = tf} < ■ ■ ■ < t™ = 1, such that 
17r n | —> 0 as n —>• oo, it therefore suffices to argue that the r.h.s. in (3.5) converges 
to the r.h.s. of (3.12). The static term in (3.12) coincides with the static term in 
(2.17). Hence, the convergence follows by use of the same arguments as in the proof 
of Theorem 2.8. Next, consider the first dynamic term in (3.5). For each cadlag path 
x e .D([0,1],M), let r m (x) := inf{7 > 0 : x* > to}. Due to the right-continuity of x, 
we have x r m( x ) > to. In consequence, 


to - Xt, 


^ ' l{xf <m<xf } - 

n—too L —' L t i -1 — TO — C/. 

i= 1 


TYl Xjmfy) 

-< 0. 

TO. — Cr m (x) 


Next, consider the second dynamic term in (3.5). First, we argue its convergence under 
assumption ii). To this end, we rewrite it as follows: 



Denote the integrand on the r.h.s. by H n . Note that the H n are predictable and uni¬ 
formly bounded. Further, since Cs = YlT=oCJP^-(t k ,t k+1 ]( s )i we ma Y choose a sequence 
of partitions ir n : 0 = < ■ ■ ■ < t™ = 1 such that H n —> H on _D([0,1],R), with 

1 {m<x t *_;Ct<xt-} m — ( t 

It follows that Jo 1 H?dM t —t f ( J HfdMt in probability; cf. Theorem 1.4.40 in [21]. 
Hence, convergence holds a.s. along a subsequence and we conclude. Under assumption 
i), the result follows by, first, integrating the pathwise inequality w.r.t. d(f>{m) and, 
then, modifying the argument along the same lines as in the proof of Theorem 2.8 (cf. 
(4.17)). □ 


4.5 Proof of Proposition 2.10 

W.l.o.g., let to. > 0 fixed and (j>(x) = 1 r x > rn \■ The general case then follows by 
integration with respect to d(j){m). The proof is based on the path-wise inequality 
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(3.5). Since is non-decreasing and M is cadlag, it implies that 


n 

1=1 


( {M U -<u) + 
V m - Cn 


(M u -Cu +1 ) + 

™ - C u+i 



n—1 

i =1 


Mt i+1 ~ M ti 
m ~ Cti+i 


n 


i=l 


m - 

m-Cu ' 


We proceed by taking expectation on both sides of this inequality and, then, passing 
to the limit. To this end, note that the expected value of the dynamic terms is bounded 
from above by zero (cf. Proposition 3.2 in [13]). Since the l.h.s. of the inequality is 
independent of the partition, it follows that 


where 


E 




< 


E [(Mi - Ci-)+] 
m - Cm 


lim [ fn(t)dQ, 

rm-ooJo 


fn(t) 


y 1 t<=(u,t i+1 ] / E [(M u ~ Cm +1 )+] E [(M ti -Cu) + ] \ 
Cti+i - Cu V m ~ Cti+i m - C u ) 


Let 


m 


E[(M t -C t ) + ]-P[M t >Ct](™-Ct) * nc 

f G c 

1 C t+ )+] E[(M t —C t )+A 

Ct+-Ct ^ m-c t + m-C* 


(4.18) 


Observe that since M t is integrable, C — > E [{M t — C) + ] /(m — C) is Lipschitz on 
(—oo ,m — e], for all e > 0. Hence, it is differentiable almost everywhere and it fol¬ 
lows that fn(t) converges point-wise to f(t), for t G [0,1] \ D, where D = fl {t : 
F(-,t ) discontinuous at Ct}, with F(x]t) := P [Mt > x\. Moreover, there is K > 0 such 
that |/ n (t)| < K , t G [0,1], n > 0. Hence, by use of dominated convergence, it follows 
that f Q f n (t)d(t —> f Q f(t)d(t and we conclude. □ 


A Appendix 

We provide a characterization of the cr-field on the canonical space fl := fl xV r . 

Lemma A.l. The Borel cr-field of the Polish space H is given by Foo := \J t>0 Ft- 
Moreover, Ft- := Vo<s<t coincides with the cr-field generated by all bounded con¬ 
tinuous functions £ : fl —> R which are Ft-measurable. 

Proof, (i) We first prove that F ^ is the Borel cr-field of the Polish space fl. Define V t + 
as the cr-field on Yf, generated by all sets of the form {9 G Yfi,9 u < s} for u G [0,1] 
and s < t; and V+ := Ui>oV t + . Then F^ = iffi <g> V+,, where F^ := Ut>oFf is the 
Borel cr-field of H (see .e.g. the discussion at the beginning of Section 1.3 of Stroock 
and Varadhan [30]). So it is enough to check that V+ is the Borel cr-field B(Yfi) of 
the Polish space Vfi. First, by the right-continuity of 0 G V+, the Levy metric on Yf 
can be defined equivalently by 

d(9,9') := inf |e > 0 : 9 t - e - e < 9' t < 9 t+e + e, Vf G Q fl [0,1]|, 
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where Q is the collection of all rational numbers. Then it follows that B(Y+) C (?{T U : 
u € Q) C V+. On the other hand, for every u € [0, 1), the map 9 H > | 9(s)ds is 
continuous under Levy metric and hence Borel measurable. Letting e —> 0, it follows 
that 9 i->- 9{u) is also Borel measurable, and hence V+, C B(Y+). We then obtain that 
Too = B(U). 

(ii) We now consider the cr-field generated by bounded continuous functions. First, it 
is well known that the filtration F° on Q is left-continuous and Ff_ = Ff is generated 
by all bounded continuous functions £i : Q —> R which are F^ continuous. 

Next, we notice that for every t > 0, 

V+ := \/ == °{TuM :« e [0,1]). 

0 <s<t 

Let £2 : V+ -> M be a bounded continuous function which is also V t + -measurable. 
Then 6((^n)«e[o,i]) = A «) r€ [o,i]), which is cr(T u At : u G [0.1])-measurable. 

On the other hand, the function 0 i— > f f“ +£ (T£(9) A t)dl from V+ to R is continuous 
and V t + -measurable for e > 0. The by taking e -A 0, it follows that T u At is measurable 
w.r.t. the cr-field generated by all bounded continuous functions £2 : Y+ -A R which 
are V f + -measurable. Therefore, is the cr-field generated by all bounded continuous 
functions on V+ which are V t + -measurable. 

Finally, since T s = <S> Vj~, it follows that Ft- = Ff_ <S> V+_. We hence conclude 

that Ft- is the cr-field generated by all bounded continuous functions £ : ft -A R which 
are ^-measurable. □ 
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